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O ! The Hartree-Fock-Bogoliubov (HFB) equations in coordinate represen- 

tation are solved exactly, i.e., with correct asymptotic boundary conditions 
£f~) ' for the continuous spectrum. The calculations are performed with effective 

Skyrme interactions. The exact HFB solutions are compared with HFB cal- 
{SJ ■ culations based on box boundary conditions and with resonant continuum 

Hartree-Fock-BCS (HF-BCS) results. The comparison is done for the neutron- 
\ rich Ni isotopes. It is shown that close to the drip line the amount of pairing 

correlations depends on how the continuum coupling is treated. On the other 
hand, the resonant continuum HF-BCS results are generally close to those of 
HFB even in neutron-rich nuclei. 
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I. INTRODUCTION 



o 

73 . 

The physics of exotic nuclei close to the drip lines has triggered a new interest for the 
study of pairing correlations in finite systems. A peculiarity of pairing correlations in weakly 
bound nuclei is their sensitivity to the effects of unbound single-particle states. 

The pairing correlations in the presence of continuum coupling have been treated both in 
HFB and HF-BCS approximations. In the HFB approximation the continuum 

is generally included in spherical systems by solving the HFB equations in coordinate rep- 
resentation. The calculations are done either in the complex energy plane by using Green 
function techniques UH, or on the real energy axis ||U . In the latter case the HFB equa- 
tions are usually solved by imposing box boundary conditions, i.e., the HFB wave functions 
are assumed to vanish beyond some distance which is chosen to be typically a few times the 
nuclear radius. 

The effect of the resonant continuum upon pairing correlations was also studied in the 
framework of the BCS approximation, both for zero l]] and finite temperature |12 . 



For deformed systems working in coordinate representation is much more difficult ||. In 
most of the deformed HFB calculations the continuum is discretized by expanding the HFB 
wave functions on a single-particle basis. Usually a harmonic oscillator basis is taken and 
one can improve the description of physical quantities at large distances like density tails by 
performing a local scaling transformation 
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The aim of this article is to show how the coordinate space HFB equations can be 
actually solved in the case of spherical symmetry and Skyrme type forces by treating the 
continuum exactly, i.e., with correct boundary conditions, and to analyse to which extent 
different approximations, namely box HFB and resonant continuum HF-BCS calculations, 
compare with the exact solutions. In this paper we will treat the continuum exactly only 
for the neutrons of neutron-rich systems for which one expects that the continuum plays an 
important role close to the drip line. For the protons we will treat the continuum by a box 
discretization. It will be shown that nuclear properties related to pairing correlations are 
correctly predicted by discretized continuum methods away from drip line but they deviate 
appreciably from exact continuum results when one approaches the drip line. 

The paper is organized as follows. In Section II we give a brief reminder of the HFB 
equations in coordinate representation. In Section III we present the procedure we have used 
for calculating the continuum HFB solutions and we discuss, in a schematic model, how the 
quasiparticle resonant states are identified. In Section IV we present the continuum HFB 
calculations for Ni isotopes in comparison with box HFB and resonant HF-BCS calculations. 
Conclusions are drawn in Section V. 



II. HFB EQUATIONS IN COORDINATE REPRESENTATION 

The HFB approximation in coordinate representation has been discussed quite exten- 
sively in the literature [@-|3J] and therefore we recall here only the basic equations. 
The HFB equations in coordinate representation read ||: 

3 (h{ra,v>a>) ~h(ra, rV) \( ^(E, rV) \ 
J \ h(ra,r'a') -h{ra,r'a') J \$ 2 (E,r'a') ) 



E + X W $i(£,ra) 
E-X l$ 2 (£,r(x) 



(1) 



where A is the chemical potential, h and h are the mean field and the pairing field, and 
represents the two-component HFB quasi-particle wave function of energy E. The mean 
field operator h is a sum of the kinetic energy T and the mean field potential T, 

h(ra, rV) = T(r, r')<W + T(ra, rV) . (2) 

The mean field potential T is expressed in terms of the particle-hole two-body interaction 
V and the particle density p in the following way: 

r(r„,rV) = /AA £ V {m ;SS,r 2 «M^,w) . (3) 

(71 <T2 

Similarly the pairing field h is expressed in terms of the pairing interaction V pa i r and the 
pairing density p: 

h(ra, rV) = / d 3 r 1 d 3 r 2 2v'v' 2 V pair (ra, r' - a'; r 1 a 1 , r 2 - cr 2 )p(r 1 a 1 , r 2 a 2 ) . (4) 

0lcr 2 

The particle and pairing densities p and p are defined by the following expressions: 
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p(ra, rV) = £ $ 2 (£„, r<r)$* 2 (E n , rV) + / dE$ 2 (E, ra)$* 2 (E, rV) , (5) 

0<E n <-\ J A 



p(r<7,rV) = £ $ 2 (£ n ,r<r)$j;(£ n ,rV) + / Cl "~ 0// dE& 2 (E, ra)^{(E, r'a') , (6) 

where the sums are over the discrete quasiparticle states with energies \E\ < —A, and the 
integrals are over the continuous quasiparticle states with energies \E\ > —A. The HFB 
solutions have the following symmetry with respect to E: 

$i(--E,ro-) = $ 2 (E,ra) () 
<f> 2 (-E,ra) = -$i(#,r<7) [> 

As it appears clearly from Eqs.([5]) and (^) we choose to work with the positive energies. 

The particle-hole and pairing interactions in Eqs.@ and @ are chosen as density- 
dependent contact interactions, so that the integro-differential HFB equations reduce to 
coupled differential equations. The zero-range character of the pairing interaction is the 
reason why one has to adopt an energy cut-off as seen in Eqs.(|5|) and (||). 

In this paper we consider systems with spherical symmetry. In this case the wave func- 
tions are readily decomposed into their radial and spin-angular parts ||: 

*i(E,T*) = Ui (Elj, r)^7(r,a) , i= 1,2 , (8) 

where: 

y%(r,a) = Y lmi (Q, $)xi/ 2 (m CT )(im z -m CT |jm) . (9) 

In the following we use for the upper and lower components of the radial wave functions the 
standard notation uij(E,r) and vij(E,r). 

As we have already mentioned, the HFB equations are usually solved by imposing to 
the radial wave functions the condition that they vanish beyond a given distance R (box 
radius). In this case the continuous spectrum is replaced by a set of discrete energies, whose 
density depends on the box radius. In what follows we discuss how the HFB equations can 
be solved by keeping the correct asymptotic conditions for the neutron wave functions. 



III. THE TREATMENT OF QUASI-PARTICLE CONTINUUM 

A. Asymptotic behaviours 

The asymptotic behaviour of the HFB wave function is determined by the physical con- 
dition that, at large distances the nuclear mean field T(r) and the pairing field A(r) vanish. 
This condition requires an effective interaction of finite range and finite-range nonlocality. 
Outside the range of mean fields the equations for &i(E,ra) are decoupled and one can 
readily find the asymptotic behaviour of the physical solutions at infinity M. Thus, for a 
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negative chemical potential A, i.e., for a bound system, there are two well separated regions 
in the quasiparticle spectrum. 

Between and —A the quasiparticle spectrum is discrete and both upper and lower 
components of the radial HFB wave function decay exponentially at infinity. For neutrons 
this implies that those components have the form, 

uij(E,r) = Ahi + \iair) , 

v lj (E,r) = Bh<i + \ip 1 r) , (10) 

where are spherical Haenkel functions, a\ = — jft(\ + E) and (3\ = — ^(X — E). These 
solutions correspond to the bound quasiparticle spectrum. In this case, the solutions are 
normalised to unity. 

For E > — A the spectrum is continuous and the solutions are: 

uij(E,r) = C[cos(5ij)ji(air) - sin(8ij)ni{a]r)\ , 

v lj (E,r) = D 1 ht\ip 1 r) , (11) 

where ji and n\ are spherical Bessel and Neumann functions repectively and is the phase 
shift corresponding to the angular momentum (Ij) . One can see that the upper component of 
the HFB wave function has the standard form of a scattering state while the lower component 
is always exponentially decaying at infinity. 

The asymptotic form of the wave function should be matched with the inner radial wave 
function, which for r — > can be written as follows: 

The HFB wave function is normalized to the Dirac 5-function of energy. This condition 
fixes the constant C to the value: 



c = tew-- (13) 

V vr n cui 

The radial wave functions are calculated by integrating the HFB equations outwards 
starting from the initial conditions (12), and inwards starting from (10) or (11) depending 
on the value of E. The solutions are propagated by a modified Numerov method towards the 
matching point, where the continuity conditions for the wave functions and their derivatives 
are imposed. These conditions determine the coefficients D ± , D 2 , D 3 and the phase shift S 
for a quasiparticle state in the continuum; in the case of a discrete quasiparticle state the 
continuity conditions and the normalisation condition determine the coefficients A, B, D 2 , 
D% and the energy E. 

The difficulty of an exact continuum calculation, i.e., with asymptotic solutions given by 
Eq.(ll), is to identify the energy regions where the localisation of the wave functions changes 
quickly with the quasiparticle energy. These are the regions of quasiparticle resonant states. 

In HFB the quasi-particle resonant states are of two types. A first type corresponds 
to the single-particle resonances of the mean field. The low-lying resonances of the mean 
field located close to the particle threshold are very important in the treatment of pairing 
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correlations of weakly bound nuclei because they become strongly populated by pairing 
correlations. 

A second kind of resonant states is specific to the HFB method and corresponds to the 
bound single-particle states which in the absence of pairing correlations have an energy 
e < 2A . In the presence of the pairing field these bound states are coupled with the 
continuum single-particle states and therefore they acquire a width. The positions and the 
widths of these HFB resonances are related to the total phase shift, calculated from the 
matching conditions, as 0: 

5(E) ~ 5 (E) + orctg , (14) 

where Er and V are the energy and the width of the resonant quasiparticle state. The 
function 5q(E) is the phase shift of the upper component of the HFB wave function in the 
HF limit, i.e., h = 0. In this limit one has: 

M>° = (E + A)$? . (15) 

If there is no single-particle resonance close to the the energy E + A in the HF limit, then 
the HF phase shift 5o has a slow variation in the quasiparticle energy region. In this case the 
derivative of the total phase shift has a Breit-Wigner form, which can be used for estimating 
the position and the width of the quasiparticle resonance. 

Thus, in the first step of the calculations we study for each channel the behaviour 
of the phase shift and we estimate the energies (widths) of the resonant states from the 
energies where the derivative of the phase shift is maximum (half of its maximum). Then, 
we choose for the integration in the energy region of a resonant state an energy grid with a 
small step. In the calculations presented below the energy step in the region of a resonance 
is T/10 and the energy cut-off is chosen to be minus the depth of the mean field. 



B. Quasiparticle resonances in a schematic model 

In order to illustrate how one can identify the quasiparticle resonances in HFB calcula- 
tions, we take here a simple model J{J0. Let us assume that the mean field is given by a 
square well potential of depth Vq and radius a. The pairing field is taken also as a constant 
inside the same radius a and zero outside. In addition, we suppose that the chemical poten- 
tial A is given. For such a system the radial HFB equations inside the potential well, i.e., 
for r < a, are: 

,1 d 2 1(1 + 1) 2 2 

\--n r 2 — + a ) u y - 7 % = o , 

J™ (XT' y£ 

A <f 1(1 + 1) „o, Or, /-.^ 

(-^r-^-^ + /3 2 H- 7 2 ^ = 0. (16) 

where a 2 = f(X + E + U ), (3 2 = ^(X-E + U ), 7 2 = f?A and U = -(V + VJ.s). 
The solutions of Eqs.(|l6D for any value of the quasiparticle energy are: 

uij = A + ji(k + r) +A_ji(k_r) , 

v tj = A + g + ji(k + r) + A^g_ji(k_r) , (17) 
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where ji are spherical Bessel functions, k± = ^t(Uq + A ± (E 2 — A 2 ) 1//2 ) and g± — (E ± 
(E 2 - A 2 ) 1/2 )/A. 

Outside the potential well the HFB equations are decoupled. In this case the type of 
solutions depends on the quasiparticle energy. They have the forms given by Eqs.(lO), (11). 

In order to simulate the potential corresponding to a heavy nucleus close to the drip line, 
we take for the model parameters the following values: Vo = 45.35 MeV, V so = 0.5 MeV, 
a = 5.2 fm, A = 1 MeV and A = -2.0 MeV. 

Here, we discuss only the quasiparticle resonant solutions induced by the bound single- 
particle states which are specific to the HFB approximation. As a typical example we take 
the case of pi/ 2 states. In the HF limit, i.e., A=0, there are two bound states at energies 
e\ = —32.873 MeV and e 2 = —10.698 MeV. When the pairing field is switched on these states 
become quasiparticle resonant states at energies E\ = 30.889 MeV and E 2 = 8.735 MeV 
with corresponding widths Fi = 0.40 keV and T 2 = 24.38 keV. These values are obtained 
by solving the HFB equations in the complex energy plane with outgoing wave boundary 
conditions. On the real energy axis one should find these two resonances from the phase 
shift behaviour. In Fig. 1 we show the phase shift (top) and its derivative (bottom) in the 
energy region of the second resonant state. One can see that the derivative of the phase 
shift is maximum at the resonance energy, and it drops to half of its maximum value when 
the energy increases by about 25 keV. This shows that the behaviour of the phase shift as 
a function of the real energy E gives accurate information on the positions and widths of 
the quasiparticle resonances. From Fig. 1 one can also see that the total phase shift does 
not cross tt/2 at the resonance energy. As discussed above, the value of the phase shift 
associated with the resonance energy is actually 8 r = tt/2 + 5q. In this case 5q = 1.59, so 
that the resonance appears when the total phase shift crosses a value close to tt and not to 
tt/2. Thus, in order to identify the resonances one can calculate the derivative of the total 
phase shift and search for the local maxima, or calculate the HF phase shift £0 and search 
for the energies associated to 5 r = tt/2 + 5 . For the 2pi/ 2 state analysed here, the two 
procedures give exactly the same position of the resonance, but this is not generally the case 
even for a square well potential fL3 |. In the present calculations we localise the resonances 
by using the derivative procedure. 



IV. RESULTS FOR NI ISOTOPES 



In this Section we apply the continuum HFB method to the case of Ni isotopes, which 
have been investigated extensively both in non-relativistic H and relativistic Hartree- 
Bogoliubov approximation [ |l4l , |l5|l . 

For the Hartree-Fock field we use the Skyrme interaction SIII whereas in the pairing 
channel we choose a density- dependent zero-range interaction: 



V = Vn 



' P{r) 



<$(ri - 1*2) , 



(18) 



with the following parameters ||: Vo = — 1128.75MeV, po = 0.134/m~ 3 and 7 = 1. 

Let us first examine the quasiparticle resonant states for the isotope 84 Ni. After conver- 
gence of the self-consistent procedure the chemical potential is A = —1.104 MeV. In Table 1 
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we show the resonant quasiparticle energies and the widths calculated from the derivatives of 
the phase shift. The quasiparticle states 2g? 3 / 2 , l<?7/2 and lh u / 2 originate from single-particle 
resonances while all the others are related to bound states. 

As already discussed in the case of the schematic model, the positions of some resonances 
may appear for values of the total phase shifts which are quite far from tc/2. We take here as 
an example the quasiparticle resonance state corresponding to the bound state 2pi/2, which 
was also analysed in the schematic model. The resonance energy and the width estimated 
from the derivative of the phase shift are E = 7.965 MeV and T = 338 keV. The value of 
the HF phase shift is in this case 5q = 0.656 so that the total phase shift associated with the 
resonance should be 5 r ~ 7r/2 + 0.656. The energy corresponding to this phase shift is E = 
7.707 MeV, which is smaller than the corresponding value extracted from the maximum of 
the derivative of the phase shift. This shows that in this case the HF phase shift has a 
non-negligible variation in the energy region of the resonance. However, in practical HFB 
calculations a small shift in the actual position of a resonance induced by the variation of 
Sq is not essential because this information is used only to fix an appropriate energy grid for 
the energy integration. 

A special behaviour can be noticed for the resonant continuum in the s\/2 channel. As 
can be seen in Fig. 2, the occupancy in this channel increases starting from —A up to an 
energy equal to 1.276 MeV. Therefore, in this channel one needs to use a very small energy 
step close to —A in order to get a correct description of the pairing correlations. Finally one 
should stress the fact that the contribution to the pairing correlations of this pronounced 
resonant structure close to the quasiparticle continuum threshold is just the manifestation of 
the loosely bound single-particle state 3si/2, which in the HFB approach is embedded in the 
continuum. This structure has nothing to do with the contribution of the S\/2 single-particle 
background continuum close to zero energy, which remains very small. 

A. Comparison between continuum and box HFB calculations 

In this subsection we analyse the sensitivity of the HFB results to the continuum treat- 
ment in the vicinity of a drip line, by comparing the results provided by continuum and box 
HFB calculations for the chain of neutron-rich Ni isotopes. The energy cut-off is the same 
in both calculations. For all box calculations presented below the box radius is taken equal 
to 22.5 fermi. Recently, some box calculations have been reported for carbon isotopes with 
box radii up to 400 fermi [TO. If these large box HFB codes could also be used for heavier 
nuclei like the Ni isotopes the differences that we show here between box and exact results 
near the drip line might be somewhat reduced. 

Lets us first discuss the properties directly related to the pairing correlations, i.e., pairing 
correlation energies and pairing densities. 

The pairing correlation energies are estimated by the difference between the total energies 
calculated in HFB and HF approach, 

E p = E(HF) - E(HFB) . (19) 

The results for continuum and box HFB calculations are shown in Fig. 3 for all Ni isotopes 
starting from A=74 up to A=88, which is the last nucleus with positive two-neutron sepa- 
ration energy, as predicted by the continuum HFB calculations (see below). Up to 86 Ni the 



7 



quantity E(HF) does not depend on the continuous single-particle spectrum. The isotope 
88 Ni is not bound in HF and therefore the E(HF) used for estimating the pairing correlation 
energy is calculated by using a box, as in box HFB calculations. From Fig. 3 one can see 
that the box HFB calculations start to overestimate the amount of pairing correlations in 
the proximity of the drip line. Thus, in box calculations the pairing energy for 84 Ni is about 
twice that of continuum HFB and it is still increasing for 86 Ni, where the continuum HFB 
calculations predict zero pairing correlation energy. 

These differences are reflected in the pairing densities, as shown in Fig. 4 for the isotopes 
84 Ni and 86 Ni. One can notice that the box calculations overestimate the pairing correlations 
in the surface region, where the localisation of the resonance wave functions with high (Ij) 
increases. Thus, in the box calculations the resonant states with high (Ij) located above 
the Fermi level are more strongly populated than the corresponding states calculated by 
using continuum HFB calculations. As an illustration we consider the occupancy of the 
single-particle resonance g 7 / 2 . In 84 Ni this resonance is located at 3.6 MeV and has a width 
of about 25 keV. If we take an energy interval 3.2 MeV< E <4 MeV around the resonance, 
we find that the total occupancy of the states in the box which are within this interval 
is about 2% higher than the corresponding occupancy in the continuum calculations. In 
box calculations the role of a resonant state is usually taken by one state with an energy 
close to the energy of the resonance, and this state has maximum localisation inside the 
nucleus. Thus, while in box calculations the pairs can virtually scatter mainly to that state 
with maximum localisation, in continuum HFB calculation the pairs can also scatter to the 
neighbouring states whose wave functions are less concentrated inside the nucleus. As a 
result the occupancy of a resonance in continuum HFB is smaller. This effect, induced by 
the width of resonant states, is missing in box HFB calculations. 

Let us consider now the two-neutron separation energies S 2n : 

S 2n = E(Z,N)-E(Z,N-2) , (20) 

which are plotted in Fig. 5. One can see that in both calculations the change of the 
sign of the two- neutron separation energies, i.e., the position of the two-neutron drip line, 
is between 88 Ni and 90 Ni, with a faster drop in the case of continuum HFB. The values 
of S 2n , evaluated within the two HFB calculations are in better agreement one with the 
other than the corresponding values of the pairing correlation energies. This is because the 
differences observed in the pairing correlation energies are much reduced when one calculates 
the differences appearing in S 2n . For the same reason one can see that even a HF calculation 
gives quite resonable values for the two- neutron separation energies close to the drip line. 
The largest differences between HFB and HF calculations appear across the doubly magic 
isotope 78 Ni. In this case the pairing energy changes quickly when two neutrons are removed 
from 1(79/2 or added to 2d 5 / 2 . Because the hole state has larger degeneracy than the particle 
state, the pairing correlations are stronger in 76 Ni than in 80 Ni. This explains the asymmetry 
seen in the behaviour of S 2n across the doubly magic nucleus 78 Ni. The fact that the value 
of S 2n predicted by HFB for 76 ~ 78 Ni is close to the data extrapolated from lighter isotopes 
indicates that the pairing interaction used in the calculations is quite reasonable, at least 
for the valence shell N=28-50. 

Next, we compare the results given by the two HFB calculations for observables related 
to mean field properties. In Fig. 6 the particle density for the isotope 86 Ni is shown. One 
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notices that the particle densities are practically the same except in the region near the box 
radius. The fact that the two particle densities are very close up to very large distances 
implies that the neutron root mean square radii (rms) calculated within the two approaches 
should be similar. This can be seen in Fig. 7 for the isotopes 80_90 Ni. 

In Fig. 7 the HF radii are also shown. In 84 Ni we can see that the HFB radius is 
slightly larger than the HF value, which is the trend usually expected when the pairing 
interaction is switched on. In this case the HFB radius is increased because the pairing 
interaction scatters some neutrons from 2d 5 / 2 to the loosely bound state 3si/ 2 which is a 
state less localised inside the nucleus. On the other hand, as seen in Fig. 7, the effect of 
pairing correlations on the radius of 86 Ni is opposite. Here, the pairing interaction scatters 
particles out of 3s±/2 state which is completely occupied in HF. The particles are scattered 
in the continuum single-particle states, mainly to single-particle resonances which have a 
larger localisation inside the nucleus than the 3si/ 2 state. Thus, in this case the radius is 
decreased when the pairing correlations are swiched on. This effect of the pairing interaction 
on nuclear radii is sometimes called "anti-halo" |6|,|T7 . 



B. Comparison between HFB and HF-BCS approximation 

The HF-BCS approximation is obtained by neglecting in the HFB equations the non- 
diagonal matrix elements of the pairing field. This means that in the HF-BCS limit one 
neglects the pairing correlations induced by the pairs formed in states which are not time- 
reversed partners. 

The extension of BCS equations for taking into account the continuum coupling was 
proposed in Refs. P,ID,|T2|. For the case of a general pairing interaction the BCS equations 
read |IJ: 



A * = E VatfijVj + £ V a ^ j g v {t)u v {e)v v {t)de , (21) 
= £ K^jfUM + £ V ue - ytvi — f f gAt'M'WVW , (22) 



^ = E^ + E/ r 9u{e)vl{e)de . (23) 

Here Aj is the gap for the bound state i and A u is the averaged gap for the resonant 
state v. The quantity g v {e) = 2 -'^ i +1 ^ is the continuum level density and 5 U is the phase 
shift of angular momentum (l u j u ). The factor g u {t) takes into account the variation of the 
localisation of scattering states in the energy region of a resonance ( i.e., the width effect) 
and becomes a delta function in the limit of a very narrow width. In these equations the 
interaction matrix elements are calculated with the scattering wave functions at resonance 
energies and normalised inside the volume where the pairing interaction is active. The BCS 
equations (21-23) are solved iteratively together with the HF equations. The corresponding 
equations are called below the resonant continuum HF-BCS equations. For more details see 
Ref. IffOl. 
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In the case of Ni isotopes the effect of the continuum is introduced through the first three 
low-lying single-particle resonances, i.e., ^3/2, #7/2 and hu/ 2 . These resonances form together 
with the bound states 2d 5 / 2 and 3si/ 2 the equivalent of the major shell N = 50 — 82. The 
energy integrals in BCS equations (21-23) are performed for each resonance in an energy 
interval defined such that |e — e u \ < 2T U , where e„ is the energy of the resonance and Y u is 
its width. In the resonant continuum HF-BCS calculations we use the same interaction as 
in HFB approach. 

In Fig. 8 we show the pairing correlation energies predicted by the resonant continuum 
HF-BCS approximation in comparison with the continuum HFB results. One can see that 
the HF-BCS results follow closely the exact HFB values up to the drip line. This shows that 
in order to estimate the pairing correlations one needs to include from the whole continuum 
only a few resonant states with their widths properly considered. 

In order to see the effect of the widths of resonant states upon pairing, we replace in 
the resonant continuum HF-BCS equations the continuum level density by delta functions. 
This means that the resonant state is replaced by a scattering state at the resonance energy, 
normalized in a volume of radius R. For this radius we take the same value as in box HFB 
calculations, i.e. R=22.5 fm. As it can be seen from Fig. 9, the pairing correlations increase 
when one neglects the widths of the resonances and the results follow closely those of box 
HFB calculations. Thus, the overestimation of pairing correlations due to the continuum 
discretisation is similar in HF-BCS and HFB calculations. 

In Fig. 7 we show also the radii calculated in the resonant continuum HF-BCS approx- 
imation. One notices that the HF-BCS radii are closer to the HF values than to the HFB 
ones. The same behaviour is found for the particle densities. This can be seen in Fig. 6 
for the case of 86 Ni, which is the last bound nucleus in the HF approximation. From Fig. 6 
one can see also that the tail of the density in resonant continuum HF-BCS calculations is 
mainly given by the particles distributed in the bound states 2<i 5 / 2 and 3«i/2 and not due 
to the particles scattered to positive energy states. In HFB calculations a part of the parti- 
cles from the bound states 2rf 5 / 2 and 3si/2 are scattered to other states, mainly to resonant 
states, with wave functions concentrated inside the nucleus. Therefore the HFB density has 
a smaller tail at large distances. 

As we have already mentioned, in the present resonant continuum HF-BCS calculations 
we neglect all the continuum contribution except for the three low- lying resonances (I3/2, #7/2 
and /in/2- This model space seems sufficient for a proper evaluation of pairing correlation 
energies up to the drip line. The rest of the continuum changes mainly the particle distribu- 
tion. In order to get a particle density closer to the HFB results one needs to introduce in 
the resonant continuum HF-BCS calculations aditional relevant pieces from the continuum. 
This work is in progress. 

V. CONCLUSIONS 

In this article we have discussed how one can actually solve the HFB equations with 
proper boundary conditions for the continuous spectrum and we have shown, for the case of 
neutron-rich Ni isotopes, how different treatments of the continuum can affect the pairing 
correlations. It was found that in the vicinity of the drip line pairing correlations are 
overestimated by the continuum discretisation done in box HFB calculations. On the other 
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hand, we have shown that the particle densities and the radii are rather insensitive to the way 
in which the continuum is treated in HFB calculations. This means that the quantities that 
are mainly related to the mean field properties do not practically depend on the different 
treatments of continuum. We have also shown that the position of the two-neutron drip 
line is not affected by the way in which continuum is treated. This is due to the fact that 
the differences observed for the pairing correlations energies in the two HFB calculations 
are diminished when the two-neutron separation energies are calculated. Moreover, the two- 
neutron separation energies predicted by HF are not very different from the HFB results. 
This shows that these quantities are not indicated for testing the pairing correlations close 
to the drip line. 

We have also analysed how the exact HFB solutions compare to the resonant continuum 
HF-BCS approximation PJTOf. It was shown that the resonant HF-BCS calculations which 
include only the first three low-lying resonances provide a very good description of pairing 
correlation energies up to the drip line. On the other hand, in the vicinity of the drip line 
the radii predicted by the resonant continuum HF-BCS calculations are larger than the HFB 
radii and closer to the HF results. This shows that one should add to the first three low-lying 
resonances additional contributions of the continuum in order to evaluate better the particle 
densities for nuclei close to the drip line. 
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TABLES 

TABLE I. Hartree-Fock single-particle energies e, HFB quasiparticle resonance energies (E) 
and widths (L) in the nucleus 84 Ni, for the various (Ij) states involved. 
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FIGURES 



FIG. 1. Phase shift (top) and its derivative (bottom) in the pi/ 2 channel for a square well model. 

FIG. 2. Occupation probability profile in the S1/2 channel for 84 Ni. 

FIG. 3. Pairing correlation energies for Ni isotopes calculated in HFB approximation. 

FIG. 4. Neutron pairing densities in in HFB calculations in 84 Ni (a) and 86 Ni (b). 

FIG. 5. Two-neutron separation energies in HFB, HF-BCS and HF approximations. For 76 Ni 
and 78 Ni the corresponding values extrapolated from experimental data [18] are also shown. 

FIG. 6. Neutron particle densities in HFB, resonant continuum HF-BCS and HF approxi- 
mations for 86 Ni. The density represented by the dotted line (HF-BCS bound) is calculated by 
including only the contribution of bound states. 

FIG. 7. Neutron rms for Ni isotopes in HFB, resonant continuum HF-BCS and HF approxi- 
mations. 

FIG. 8. Pairing correlation energies calculated in resonant continuum HF-BCS approximation 
compared to continuum HFB. 

FIG. 9. Pairing correlation energies calculated in the resonant continuum HF-BCS approxima- 
tion by neglecting the widths effect compared to box HFB results. 
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Occupation probability in the s1/2 channel 
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